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The superconducting state of an optimally doped single crystal of Ba(Feo.93Coo.o7)2As2 was in- 
vestigated by '^^As NMR in high magnetic fields from 6.4 T to 28 T. It was found that the Knight 
shift is least affected by vortex supercurrents in high magnetic fields, > 11 T, revealing 

slow, possibly higher order than linear, increase with temperature at T < 0.5 Tc, with Tc ~ 23 A". 
This is consistent with the extended s-wave state with Aig symmetry but the precise details of 
the gap structure are harder to resolve. Measurements of the NMR spin-spin relaxation time, 
T2, indicate a strong indirect exchange interaction at all temperatures. Below the superconducting 
transition temperature abrupt changes in vortex dynamics lead to an anomalous dip in T2 consistent 
with vortex freezing from which we obtain the vortex phase diagram up to H = 28 T. 

PACS numbers: 



I. INTRODUCTION 

As a result of the discovery of pnictide superconduc- 
tors by Kamihara et al.^ in 2008 there has been intense 
interest in this new family of materials. In the past two 
years their quality has greatly improved and large sin- 
gle crystals have become available. However, the nature 
of the superconducting state and the corresponding gap 
structure of the pnictide superconductors is not settled. 

While the pairing mechanism, related to magnetic fluc- 
tuations, and basic structure of the extended s-wave 
state are generally agreed upon^"— details of possible gap 
anisotropy are currently actively investigated, often with 
opposite conclusions. 

Experiments from angle resolved photoemission spec- 
troscopy (ARPES)r'"— Andreev- reflection spectroscopy;^ 
and specific healjiS measurements support the fully- 
gapped model. But penetration depth measure- 
mentsr^"— nuclear quadrupole resonance (NQR), and 
specific heat measurementsi^ have been interpreted in 
terms of nodal gap structure. These measurements in- 
dicate that the superconducting order parameter might 
be fundamentally different in different classes of these 
materials, further complicating interpretation of the ex- 
periments. 

Moreover, even in a single compound, but with differ- 
ent amount of doping, various gap signatures have been 
seen, from fully gapped, to nodal character. For example, 
calorimetricj^"— transport^ and optical conductivit}*^ 
measurements in overdoped Ba(Fei_2:Coa:)2As2 {x > 
0.1) show nodal character, and similarly in underdoped 
samples^^i^i whereas at optimal doping (a; ^ 0.06 — 0.08) 
results are contradictorily interpreted either in terms of a 
more or less isotropic gap ^^i^^'^^ or a strongly anisotropic 
gap^ 

The fact that the gap structure might be doping- 
dependent is a possibility. Theory predicts that in multi- 



band materials, such as pnictides, one may naturally find 
states with "accidental" nodes, even in the most symmet- 
ric Aig configuration. In most cases the nodes appear on 
the outer electronic Fermi surface sheets, depending on 
the values of the interaction parameters, which are func- 
tions of doping/^"— The outstanding question now is 
what is the form and position of nodal lines on the elec- 
tron Fermi surface<22^"— Other nodal locations (e.g. on 
central hole pocket) have also been considered 1^1^ 

In this work we report nuclear magnetic reso- 
nance (NMR) measurements on single crystals of 
Ba(Feo.93Coo.o7)2As2. NMR has advantages in that it 
can probe bulk characteristics, is sensitive to the elec- 
tronic structure within a penetration depth of the sur- 
face of a single crystal, > 100 nm, and is less sus- 
ceptible to surface conditions as might be the case in 
STM and ARPES experiments. In principle. Knight 
shift data from NMR can be an appropriate indicator 
of fully gapped or nodal gap structure. For the former, 
the Knight shift has the temperature dependence of the 
Yosida function^ or normal fluid density, and for a gap 
with line nodes the temperature dependence should be 
linear at low temperaturesi^ Our goal was to investi- 
gate the Knight shift at sufficiently high magnetic fields 
that vortex contributions to the temperature dependence 
of the local field are minimized, and at sufficiently low 
temperatures to identify the gap structure. It is also 
important to determine independently that the vortex 
structures have a static distribution, that is to say there 
is a solid vortex state. Then the vortex configuration is 
stable and cannot introduce temperature dependence to 
the lineshape, affecting the Knight shift analysis. Conse- 
quently we have measured the spin-spin relaxation time, 
T2, as a function of temperature and magnetic field and 
found it to be a good indicator for the onset of vortex dy- 
namics. We focus on measurements of the Knight shift 
below the temperature where vortices become frozen. 
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II. EXPERIMENTAL METHODS 

We have performed ^^As NMR studies on a single crys- 
tal of Ba(Fei_a;CO:c)2As2 {x = 0.074) with zero field 
Tc = 22.5 K, varying temperatures from 2 K to 200 K 
with external magnetic field from H = 6.4 T to = 28 
T parallel to the c-axis of the sample. The measure- 
ments were performed at Northwestern University and 
the National High Magnetic Field Laboratory in Talla- 
hassee, Florida. The crystal was grown at Ames Lab- 
oratory by the self-flux method^ having dimensions of 
3.9x5.3x0.7 mm^ and mass of 49.3 mg. A Hahn echo 
sequence (7r/2 — r — tt — t— echo) was used for spectra, 
Knight shift, linewidth, and some of our spin-spin relax- 
ation time, T2, measurements, where r is the delay time. 
A typical tt/2 pulse length was 4 fj,s, defined as that which 
gave maximum echo intensity. We used frequency sweeps 
when a single pulse did not cover a sufficiently wide fre- 
quency range. Delay times for acquisition of the spectra 
were varied from 70 to 200 fis depending on temperature. 
The Knight shift and linewidth measurements were deter- 
mined from gaussian fits to the line shape. The gyromag- 
netic ratio of the bare nucleus ^^As, 7 = 7.2919 MHz/T, 
was used as the reference for the Knight shift. Linewidths 
were defined as full-width-at-half-maximum (FWHM) of 
the gaussian distribution. For T2 measurements from the 
Hahn echo sequence, we varied t from 100 fis to 1.40 ms, 
and determined the rate from the initial portion of the 
recovery. In high magnetic fields, H > 17 T in a re- 
sistive magnet, we used the Carr-Purcell-Meiboom-Gill 
(CPMG) sequence which is resistant to magnetic field 
fluctuations {tTx/^ ~ t — iTy — r— echo— r — tTj, — ... and 
x,y indicate orthogonal RF phases). Typically, the refo- 
cusing time, 2r, was chosen to be 250 ^s. In the course 
of these spin-spin relaxation measurements our compar- 
ison of CPMG and Hahn echo experiments revealed the 
existence of strong field fluctuations intrinsic to the sam- 
ple. In all cases, the pulse sequence repetition time in 
the superconducting state was taken to be of order the 
spin-lattice relaxation timcj^S Ti, and was increased with 
decreasing temperature. 



III. NMR SPECTRUM AND KNIGHT SHIFT 

Our fleld swept spectra of ''^As NMR and the shift 
of their peak frequencies with temperature are shown 
in Fig. [TJ From the satellite transition we find the 
quadrupolar coupling, i/q = 2.58 MHz (0.35 T), similar 
to the value reported by Ning et al^ The quadrupolar 
satellite transitions are suppressed because of disorder in 
the electric field gradient (EFG) tensor, most likely from 
Co doping, from which the central transition is largely 
immune. In part a) of this figure there is a small com- 
ponent of the spectrum having a negligible Knight shift 
(0.003%). It is clear that these nuclei are from impurity 
or defect sites having a different orbital electronic shift. 
The spectra in the lower panel are obtained from the cen- 



tral transition. Similar observations were made by Ning 
et alM 

Our measurements of the temperature variation of the 
total Knight shift, K, of the ^^As central transition are 
presented in Fig. [2] Spin-singlet pairing is evident from 
the sharp decrease of K{T) at Tc- The onset of this de- 
crease at 18 K in 11.6 T coincides with the measured 
from the onset of the drop of the resistance reported by Ni 
et alM- in the same magnetic field. This is in contrast to 
the cuprates, especially the underdoped cuprates, where 
the Knight shift starts to decrease well above Tc, indicat- 
ing the presence of a pseudogap. At low temperatures in 
the superconducting state, we observe that the Knight 
shift decreases, almost linearly with temperature. 

We express the total Knight shift as K = Ks+ Koih + 
Kq -\- Kpara- Thc Spin part of the shift is Ks = ^hfXsj 
where Xs is the spin susceptibility, and Ahf is the hyper- 
fine field. Both the orbital, i^orb, and the quadrupolar 
shift. Kg, are expected to be temperature independent. 
Additional magnetic shifts which we call Kpara could 
be associated with impurities;^ and in principle, could 
be temperature dependent, although this is not seen in 
cuprateS )^^'^^ nor do we find any evidence for this in the 
present work. 

In the normal state of Ba(Feo.93Coo.o7)2As2 below 170 
K, we find the Knight shift can be described by an acti- 
vated process. However, there is a substantial range at 
low temperatures where the Knight shift is essentially 
temperature independent. Fig. [51 We use an Arrhe- 
nius form, K{T) = A + i? exp(— Ag/T), where we find 
A, = 394 ± 16 K, A 0.252%, and B = 0.179 with 
the fit shown by the (red) curve in the figure. At high 
temperatures, 90 - 170 K, the total Knight shift is pro- 
portional to the bulk susceptibilityj^i giving a hyperfine 
field of ^HF = 14 kOe/fiB, comparable to the value of 
18.8 kOe//XB reported by Kitagawa et al^ Therefore the 
thermally activated process modifies thc local density of 
states and directly affects the spin susceptibility. How- 
ever, neither the origin of this high-temperature behav- 
ior, nor its relation to superconductivity is known. The 
first observations of this effect, along with a similar anal- 
ysis, was performed by Ning et al^ giving As = 520 K 
for their 8 % Co crystal of Ba(Feo.92Coo.o8)2As2. 

In thc superconducting state, we found that the total 
Knight shift gradually decreases on cooling from 12 K 
to 2 K in 11.6 T, figures 2 and 3. Similar behavior is 
consistently observed for higher magnetic fields up to 20 
T, Fig. 4. 

If one takes the view that this is a linear temperature 
dependence, then such behavior can be associated with 
gap nodes at the Fermi surface, similar to cuprate su- 
perconductors with c?-wave symmetryj ^^'^^i^^ The Knight 
shift Kg (T) is proportional to the susceptibility 

f,ENiE)l^^^iJ^in^I^^±J^ (1) 

Xn J 2flH 

where Xn = '2fJ,'^Nf - normal state susceptibility, f{E) 
is the Fermi-Dirac distribution function, and N{E) is 
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the density of states. For temperatures T/Tc > ijlH/2Tc 
[f[E - fxH) - f{E + fiH)]/[2fiH] « -df{E)/dE and for 
nodal quasiparticles N{E) ~ E/Aq, resulting in a lin- 
ear temperature dependence x/Xn ~ T/Tc x const 0{1). 
This is a much steeper rise than what we observe. Thus, 
we can exclude the possibility of nodes simultaneously 
present on both hole and electron sheets. However, nodes 
might be present on one of the pockets, presumably elec- 
tronic, which in combination with the fully gapped hole 
pocket would give a slower T-dependence. Another pos- 
sibility is that the observed temperature dependence is 
higher power than linear, and the electronic supercon- 
ducting gap in Ba(Feo.93Coo.o7)2As2 is either anisotropic 
with weak "accidental" nodes, or has only a minimum 
with typical energy scale Amin < fJ^H, or even isotrop- 
ically gapped, with strong pairbreaking due to impuri- 
ties.—"^ In this regard, the existence of impurities affect- 
ing the local magnetic fields is evident in the temperature 
dependence of the NMR linewidth in the normal state. 
Fig. 5. 

Comparison of our Knight shift measurements for H = 
11.6 T with the above model and different gap struc- 
tures are shown in Fig. |3l The spin susceptibility is 
calculated according to Eq. [1] where the shift of the 
energy levels in magnetic field is given by parameter 
Z = fisH/Tco ^ 0.3 {H = 11 T), where hb is the 
Bohr magneton, and Tco = 25 K is an estimate of the 
zero-field transition temperature in the absense of impu- 
rities. The density of states is computed in the two-band 
model used previously and described in references 30 and 
46. The gap on the hole FS is isotropic, Ah{(j)) = Ai, 
while on the electron FS it is a function of the angle 
Ae((/)) = A2(l — r + rcos2(/)), with r = being the 
isotropic gap, and r = 1 gives four equally spaced nodes. 
The impurity scattering is given by the concentration 
parameter 7 = (nimp/T^Nf)/2TTTco, interband potential 
fraction {SV = ^12/^11), and an intrahand scattering 
phase shift (tan (5 — TrNfVii). The zero-temperature 
shift Ks {T = 0) is not well-defined experimentally and is 
left as a free parameter. Below 10 T the temperature de- 
pendence is stronger and it is plausible that this can be 
associated with diamagnetic screening currents or vor- 
tex supercurrents which do not cancel as effectively in 
the intervortcx region at low fields. At high tempera- 
ture (0.5 < T/Tc < 0.8); the discrepancy between our 
data and the model is likely due to thermal fluctuations 
of vortices that introduce a temperature dependence to 
the NMR lineshape in a region where vortices are be- 
ing pinned. It is just in this region of temperature that 
slow fleld fluctuations from vortices are manifest in the 
linewidth. Fig. 5, and in spin-spin relaxation, discussed in 
the next section. Fig. 6. We observe a systematic trend in 
the data of increasing Knight shift, K{0), with increasing 
fleld. Fig. 4, which might be accounted for by considering 
the vortex core contributions to the density of states. 

From Fig. 3, in the temperature range T/Tc < 0.5 
where we can compare the Knight shift data with the 
theoretical model, we find that more or less all electronic 



gap structures that we have considered in the Aig sym- 
metry class can fit the data. One might say that the in- 
termediate cases of anisotropic gaps with either minima 
r = 0.45 or close nodes r = 0.55 are less suitable, espe- 
cially in the strong scattering case, leaving surprisingly a 
possibility of either an isotropically gapped electronic FS 
or one having a gap ^ cos 2(f). A more reliable identifica- 
tion of the angular structure could be possible in lower 
magnetic fields {Z = inset in Fig. 3(a)), but the analysis 
would be quite complicated due to the effects of vortices 
and diamagnetic screening currents. Different anisotropy 
parameters r in this model give a reasonable coverage of 
the phase space of possibilities for the spectrum of quasi- 
particles. We checked that with a more realistic three di- 
mensional Fermi surface and the gap suggested in Ref . [13 
results do not change much (dashed line in the inset of 
Fig. 3(a) that appear exactly between the isotropically 
gapped r — and purely nodal r = 1 cases). 

The linewidth of the central transition. Fig. [5l full- 
width-at- half-maximum (FWHM), was determined by 
fitting the spectra to a gaussian function. The flt devi- 
ates slightly near the wings but at low temperatures the 
lineshape appears to be symmetric. The linewidth in- 
creases with decreasing temperature in the normal state. 
Below Tc, there is a decrease, followed by an upturn at 
~ 11 K in 11.6 T. These features closely parallel all of 
the aspects of the NMR linewidth previously reported 
by Chen et alMi^ for ^'^O NMR on high quality sin- 
gle crystals of Bi2Sr2CaCu208+5 (BSCCO). For cuprate 
supcrconductors^ii^ this similar behavior in the normal 
state was identifled with impurities. Its classic signature 
is a Curie or Curie- Weiss temperature dependence to the 
linewidth and a constant Knight shift. The basic argu- 
ment is that impurities introduce strong polarization of 
the nearby conduction electrons which produce a spa- 
tially oscillating spin density that couples through the 
hyperflne interaction to the nuclei (RKKY- interaction). 
Consequently, the average local fleld is not signiflcantly 
perturbed, so the Knight shift is unaffected; but, a large 
distribution of Knight shifts gives a signiflcant broaden- 
ing of the NMR spectrum. The polarization of paramag- 
netic impurities is typically quite temperature dependent 
(Curie- like), and is the likely origin for the temperature 
dependence in the linewidth below T < 80 K. At higher 
temperatures there is a possible additional temperature 
dependence from the activated process that is evident in 
the Knight shift shown in Fig. 2. 

The success we have had with the above model in 
understanding the ^^O NMR linewidth in BSCCO,i^ 
serves as a guide for our interpretation of the im- 
purity contributions to the ''^As NMR linewidth in 
Ba(Feo.93Coo.o7)2As2. Wc focus on the range of tem- 
perature where the Knight shift is mainly temperature 
independent, < T < 80 K. We fit the hncwidth to the 
following phenomenological relation. 



Ai/ = Ai^o (1 - aexpi-AJkBT)) + —^K,{T). (2) 
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Since there may also be a distribution of Knight shifts 
with a thermahy-activated origin, we aUow for this contri- 
bution with the additional term, aexp(— A^/fesr), which 
goes beyond the impurity model used for BSCC04^ 
Here, = 394 K, obtained from a fit to K{T) in 
the normal state with Korb = 0.229% with parameters, 
Ai/QjC, and a. Overall, the fit is reasonable, capturing 
the essential characteristics of the normal state behav- 
ior of the linewidth with a reduction in the supercon- 
ducting state as the impurity local field distribution be- 
comes gapped following the temperature dependence of 
the Knight shift. Although impurity effects can account 
for the temperature dependent linewidth, we cannot com- 
ment on what type of impurity might be responsible. 

The upturn of the linewidth at 11 K in 11.6 

T in the superconducting state is the mark of vortex 
freezing. It occurs at T^/Tc 0.61 in fl" = 11.6 T, 
and indicates that there is a substantial temperature 
range below Tc where fluctuating vortices are in a liquid- 
like phase. We can compare this with YBa2Cu307_j, 
(YBCO)i^ and overdopcd BSCCO^i where vortex melt- 
ing for the same magnetic field occurs at T/Tc = 0.87 
and 0.21 respectively. Calculations of the melting tran- 
sition in high magnetic fields for pnictidc superconduc- 
tors have been performed by Murray and Tesanovio^. 
However, a true, thermodynamic, vortex-liquid phase in 
Ba(Feo.93Coo.o7)2As2 would be surprising since both the 
mass anisotropy^ and Tc are small ; for example, as 
compared to cupratc superconductors. It was argued by 
Yamamoto et al^ that a confluence of weak thermal fluc- 
tuations and a very broad distribution of pinning forces 
gives rise to an irreversibility field significantly less than 
the upper critical field in crystals of Ba(Feo.9Coo.i)2As2. 
Generally the NMR spectrum is sensitive to field fluctu- 
ations on time scales slower than the Larmor precession 
period which is tens of nanoseconds in our case. Con- 
sequently, the sharp variation in the linewidth that we 
observe at Tm suggests that vortex dynamics are slower 
than this for T < Tm- To investigate this further, we have 
performed spin-spin relaxation experiments discussed in 
the next section. 

Usually, the vortex contribution to the NMR or fiSR 
linewidth^i^ is interpreted in terms of the penetration 
depth. In Ba(Feo.93Coo.o7)2As2 the linewidth broaden- 
ing from the vortex solid state can be seen in Fig. [5] as a 
difference between data (represented by a black dashed 
line) and our model for the vortex free linewidth (solid 
blue curve) discussed above. The vortex component of 
our linewidth at the lowest temperatures, 7.5 kHz 
(~ 10 gauss) is much less than predicted by Ginzburg- 
Landau theory^ which is 33 kHz for a penetration depth 
from ^SR of Xab = 217 nm,— or 15 kHz for the value, 
325 ± 50 nm, from scanning SQUID and magnetic force 
measurements)^ The most likely reason for this discrep- 
ancy is c-axis vortex disorder discussed by Brandt.— 
For example, at very low magnetic fields, H < I T, 
this plane-to-planc disorder is well-established in BSCCO 
from /iSR and NMR mcasurementsiii^"— and can be 



expected in Ba(Feo.93Coo.o7)2As2 since in-plane disorder 
and strong pinning has already been established by scan- 
ning probe measurements^^^ and small angle neutron 
scattering 1^ 



IV. SPIN-SPIN RELAXATION AND VORTEX 
DYNAMICS 

To explore possible effects of vortex dynamics, already 
indicated in our NMR linewidth measurements, we have 
investigated spin-spin relaxation. There are two classes 
of spin relaxation experiments that we perform. The 
results measured by Hahn echo methods arc presented 
in Fig. 6 and 7, and by the CPMG pulse sequence in 
Fig. 6 and 8. The nuclear spin dephasing times, T2H and 
T2CPAIG-: from these experiments can be viewed semi- 
classically as a measure of field fluctuations along the 
direction of the applied magnetic field. The dephasing 
rate is a summation of contributions from spin-lattice re- 
laxation (Redfield contribution) , local fiuctuations in the 
external magnetic field, and fiuctuations in the internal 
fields such as those from the nuclear dipole-dipole inter- 
action, or from vortex dynamics. The Hahn echo method 
is generally useful for homogeneous relaxation (all nuclei 
are equivalent) in a steady external field. However, if 
the local field fluctuates inhomogeneously, or there is nu- 
clear spin diffusion in a spatially inhomogeneous steady 
field^i^ then the CPMG method can be helpful in iden- 
tifying the time scale of the fluctuations through varia- 
tion of the spin refocusing time, r (see the experimental 
methods section) 1^ 

From our relaxation measurement data on 
Ba(Feo. 93000.07)2^82 shown in Fig. 6, it is immedi- 
ately apparent that T2CPMG is signiflcantly longer than 
T2H and that in the normal state, both are much longer 
than the temperature independent dipolar contribution, 
0.7 ms, calculated from the rigid lattice limit, shown as 
a red dashed line. For all CPMG experiments shown 
here 2t was set to be relatively short, 250 /iS to 700 /iS. 
We infer that significant dephasing takes place on time 
scales longer than 700 /iS and, since this is in the normal 
state, the process responsible has nothing to do with 
vortices. Rather it is an indication of the presence 
of slow field fluctuations from an unidentified source. 
In fact, by varying the refocusing time, 2r, we have 
determined that this dephasing is effectively eliminated 
in the CPMG results which we report here. We defer 
our discussion of the mechanism for the Hahn echo 
relaxation to later in this section. 

The temperature dependence of T2CPMG in the normal 
state can be understood from the temperature depen- 
dence of spin-lattice relaxation, through the Redfield con- 
tribution!^ This behavior can be expressed phenomeno- 
logically as, T2^Vmg = T^^ + [aTi)-^^ Taking Ti from 
Ning et al^ we find a = 5.7 and T2 = 12 ms, where T2 
is the spin-spin relaxation time after correction for this 
contribution. Since Ti is long enough below Tc (> 150 
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ms) we may take T2CPMG to be equal to T2 in the su- 
perconducting state. Nonetheless, we find that there is a 
remarkable dip in both T2CPMG a-nd in T2H at T,„, Fig. 6 
to 8, which we identify with thermal fluctuations of vor- 
tices. This follows from the fact that the temperature 
of the minimum in T2 coincides with the vortex freezing 
temperature apparent in the NMR linewidth, presented 
in Fig. 5 for 11.6 T and that this dip at appears just 
below the upper critical field, Fig. 9. 

Vortex dynamics contribute similarly in these two 
types of spin-spin relaxation experiments although the 
contribution to the Hahn echo can be an order of mag- 
nitude larger than for the CPMG, especially at higher 
magnetic fields. This means that the principal fluctua- 
tion time scale must be of order T2H: but longer than 
the refocusing time, 2r = 700 ^s, in the CPMG mea- 
surement. 

In previous work on ^^O NMR in YBCO it was found 
that the spin-spin relaxation rate, I/T2, had a peak below 
Tc^^— Kramer et al^ proposed a dynamical charge- 
density-wave state in YBa2Cu307_5 coupled to the nu- 
cleus through the quadrupolar interaction. However, 
Bachman et al^ found that this peak was an artifice of 
vortex dynamics which become slow enough in the vor- 
tex solid to be observable in T2 ■ They found a lorcntzian 
component to the spectral density for field fluctuations 
that abruptly onsets at the irreversibility temperature. 
The fact that the onset temperature was magnetic field 
dependent indicated that the phenomenon was not re- 
lated to a charge density wave. This argument holds 
equally well in the present case for Ba(Fco.93Coo.o7)2As2. 
Here, as for the cuprates, we associate the temperature 
of the maximum spin-spin relaxation rate (the dip in the 
relaxation time) with the temperature where the vortex 
dynamics rapidly slow down, i.e. the irreversibility tem- 
perature. For comparison, in Fig. 9 we plot, along with 
our data, the irreversibility curve from Prozorov et al^ 
obtained at lower magnetic fields for similar crystals from 
the Ames Laboratory. Our data includes the tempera- 
tures of anomalies from the NMR linewidth (spectrum) 
and T2 measurements we have discussed previously. The 
CPMG sequence was used to obtain T2 at the higher 
magnetic fields, 20 T to 28 T because the Hahn echo 
experiment is adversely affected by magnetic field fluctu- 
ations endemic to high field resistive magnets.— We also 
compare our measurements of the superconducting tran- 
sition with those of Ni et alM^ providing a consistent 
vortex phase diagram for Ba(Fco. 93000.07)2^82. 

Now we return to a discussion of the large T2 we 
have observed, as compared with the expected tempera- 
ture independent upper bound established by the nuclear 
dipole-dipole interaction, the red dashed line in Fig. 6. In 
order for T2 to be longer than the dipolc limit there must 
be some form of coherent averaging which in liquids is the 
well-known phenomenon of motional averaging. In met- 
als there is a similar mechanism, a type of spin-motional 
averaging, which is attributable to a strong indirect inter- 
action between adjacent nuclei, ^^As in our casci^^iii For 



this process to work the coupled nuclei must be resonant, 
that is to say they must have the same Larmor frequency 
to within the small band given by the bare nuclear dipole- 
dipole interaction. This requirement is imposed by the 
necessity to conserve energy in the spin 'flip-flop' process 
for spin-spin relaxation. The indirect interaction is a cou- 
pling between nuclei, mediated by the conduction elec- 
trons, involving an RKKY type polarization and the hy- 
perfine interaction. For metallic elements with only one 
spin isotope this leads to averaging of the dipolar local 
fields, called exchange narrowing, and can substantially 
increase T2, a factor of 5.6 in the case of platinum<^i2S 
For elements with more than one spin isotope the nuclei 
are non-resonant and the interaction leads to T2 shorter 
than the dipolar limit, as is the case for thallium. From 
the theory of exchange narrowing^ we have. 



1/I2 ^ W^poIeM (3) 

where uje is the frequency of the exchange interaction 
and LUdipoie is the bare dipolar interaction between nuclei. 
From this relation we find, uje — 24 kHz, which is larger 
than the 4 kHz,— observed for metallic platinum, but not 
unreasonably large given that ^^As is 100 % naturally 
abundant as compared to 33.8 % for ^^^Pt. Although 
this explanation can account for the longer T2 than ex- 
pected from the dipolar interaction, it does not explain 
the shorter phase coherence time found from Hahn echo 
measurements as compared with the CPMG sequence. 

There are two possibilities for the source of low fre- 
quency field fluctuations that could give rise to the dif- 
ference between Hahn and CPMG results. In the first, 
we consider fluctuating magnetic fields, either from impu- 
rities or residual effects of the suppressed antiferromag- 
netism inherent to the BaFe2As2 system. However, it is 
hard to reconcile this explanation with our observation 
that T2H is independent of magnetic field and tempera- 
ture, as is evident in Fig. 7 after allowing for the vortex 
contribution. One might expect that thermally driven, 
magnetic field fluctuations would be both temperature 
and magnetic fleld dependent over the wide range we 
have covered. A second possibility is that the rather large 
indirect interaction leads to nuclear spin diffusion that 
produces dephasing in a locally inhomogeneous, static, 
magnetic field. We have estimated the diffusion coeffi- 
cient from the indirect interaction^i^^ to be Ds ~ 10~^^ 
cm^/s which, in order to explain our Hahn echo results, 
would require an average magnetic field gradient of « 1.5 
X 10^ G/cm. Associating a magnetic field inhomogene- 
ity of this size with residual magnetism in the sample 
could provide the temperature independent mechanism 
that would be required to account for the Hahn echo re- 
sults. The length scale for the distribution of such mag- 
netism is constrained by the linewidth to be less than 
w 20 nm. Although it may seem that the latter explana- 
tion is the more likely, further work will be necessary to 
confirm or disprove cither of these suggestions. 
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V. CONCLUSION 

We have studied ^^As NMR spectra and spin-spin 
relaxation on a single crystal of optimally doped 
Ba(Feo. 93000.07)2^82 from room temperature to well 
below the superconducting transition temperature in 
magnetic fields as high as 28 T. The temperature 
dependence of the Knight shift is compatible with a 
general, two-band Aig symmetry scenario for the super- 
conducting states. The details, nodes or full gap, cannot 
be determined from the present work. Comparison with 
an elementary 2D model might indicate that the order 
parameter in this compound has nodes but they have 
relatively weak weight compared to the fully gapped 
portions of the Fermi surface. We have identified a 
sharp signature of irreversibility in the temperature 
dependence of our measured spin-spin relaxation from 
which we have established a vortex phase diagram up to 
high magnetic fields. We have found that there is a very 



strong indirect exchange interaction in this compound 
providing a possible explanation for our observations of 
spin-spin relaxation attributed to nuclear spin diffusion 
in static magnetic field gradients. 
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FIG. 2: Knight shift of Ba(Feo.93Coo.o7)2As2. The temper- 
ature dependence of the Knight shift consists of three parts: 
at low temperature (below the vortex freezing temperature), 
near the superconducting transition (the arrow indicates Tc), 
and in the normal state. The drop in Knight shift at the su- 
perconducting transition is evidence for a spin-singlet state. 
In the normal state the Knight shift follows an activated be- 
havior, with activation energy of 394 K. The solid curve (red) 
is the fit to this activated process. The low temperature be- 
havior shows a decreasing Knight shift from which we find 
K{Q) = 0.231% by extrapolating to zero temperature. 



FIG. 1: ■^'^AsNMR spectra of Ba(Feo.g3Coo.o7)2As2. a) Single 
crystal spectra were obtained from a field sweep at 84.88 MHz. 
From the quadrupolar satellite transitions, shown magnified 
for clarity, uq was found to be 2.58 MHz (0.35 T). b) Spectra 
of the central transition are shown at different temperatures, 
normalized in area. 
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T/T, T/T, 

FIG. 3: Normalized susceptibility of Ba(Feo.93Coo.o7)2 As2. 
The theoretical curves (solid lines) are shown for several pos- 
sible gap structures on the electronic pocket, given by the 
form Ae((^) = A2(l — r + r cos 20), for intermediate (panel 
a) and strong scattering (panel b). To compare these results 
with the experiment we took Z = ^bH/Tco = 0.3 [H — 11 
T). Inset in panel (a) shows results for the clean case and at 
low fields {Z — 0) - note that aside from the low-T region all 
anisotropic gaps, r — 0.45, 0.55, 1, and a model with more 
realistic 3D Fermi surface, behave in approximately the same 
way. In high field and for weak impurity scattering (a), all 
gap models are consistent with experiment at low T. In the 
unitary limit (b), the best fit is given by either combination 
of two isotropic gaps Ah{(l>) ~ —Ae{<j)) = A (r = 0), or, sur- 
prisingly, by an isotropic hole Ah(0) = Ai and d-wave like 
electron gaps Ae(</>) — A2Cos2</i (r = 1). It is possible that 
the dip in the data, 0.5 < T/Tc < 0.8, is a consequence of 
changes in the NMR lineshape owing to vortex , section IV. 
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FIG. 4: Knight shifts at various magnetic fields. The low field 
values have a strong temperature dependence at low temper- 
atures. This is ascribed to vortex supercurrent contributions 
to the line shape which are expected to become negligible at 
high fields. The weak temperature dependence of the Knight 
shifts for J/ > 11 T at the lowest temperatures, should be 
independent of vortex supercurrents. An increase in Knight 
shift, K{0), with increasing magnetic field is observed. 
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FIG. 5: Linewidth of the ^^As central transition at H — 
11.6 T. The upturn of the hnewidth at Tm ~ 11 K, indicates 
freezing of the vortex liquid^ (the inset shows an expanded 
view of this transition). The solid blue curve, down to Tm, 
is a fit to the model described in the text. Below Tm the 
solid blue curve is a calculation from this model, absent any 
vortex contributions to the linewidth. The dashed black line 
is a guide-to-the-eye representing the data. The difference at 
low temperatures between the solid and dashed curves in the 
inset is the vortex contribution to the linewidth which we find 
to be « 7.5 kHz. 



10 
E, 4 

CM 

I- 2 
1 



2 

_^ (a) Ba(Feog3COoo7)2As2 
„ _ H = 6.4 T // c-axis 



f A Hahn 

-*■ n CPMG 



50 100 150 200 250 300 



FIG. 6: Spin-spin relaxation times at H — 6.4 T. Mea- 
surements of spin-spin relaxation times with Hahn echo and 
CPMG methods are shown for comparison. A superconduct- 
ing magnet was used to rule out T2 reduction due to field 
fiuctuation. a) The slow increase with decreasing tempera- 
ture down to Tc corresponds to reduction in the Redfield con- 
tribution, i.e. from spin-lattice relaxation, b) An expanded 
view at low temperature below Tc, shows that there is a de- 
crease in both relaxation times. We associate this with vortex 
dynamics providing a maximum contribution to the rates at 
Tm, taken from Fig. 5. The red dashed line indicates T2 in 
the dipolar limit, 0.7 ms. 
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FIG. 7: Hahn echo relaxation times, T2H, for various magnetic 
fields from a superconducting magnet. There is a clearly de- 
fined temperature, Tm, for a minimum in T2 at each magnetic 
field which we associate with vortex freezing. The minimum 
T2 becomes even smaller than the T2 from the dipolar limit 
(0.7 ms) when the magnetic field increases. Additionally, T2H 
shows a decrease in the normal state just above Tc ai H = 16.8 
T. 
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FIG. 8: CPMG relaxation times, T2cpmg, for various mag- 
netic fields. The transition temperature^^ (black upward ar- 
row) and the vortex freezing temperature (red downward ar- 
row) correspond to abrupt changes in relaxation. The dashed 
black arrows indicate the minimum relaxation time from the 
Hahn echo method. 
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FIG. 9: Vortex phase diagram. The vortex freezing temper- 
ature deduced from the minimum in the spin-spin relaxation 
time and from the minimum in the linewidth is shown as a 
function of magnetic field along with the vortex irreversibil- 
ity temperature from Prozorov et al-^ The superconduct- 
ing transition temperature inferred from the NMR spectrum 
linewidth and Knight shift are compared with the report from 



Ni et al,^ 



